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I. INTRODUCTION 

Density matrix renormalization group (DMRG) 
method has been a powerful numerical tool since White's 
pioneering works Q. DMRG is a real space numeri- 
cal renormalization method by using the reduced den- 
sity matrix of a large system to approximate the ground 
(or excited) state(s) of the system with great accu- 
racy. There are many DMRG applications for (quasi)one- 
dimensional (ID) systems, e.g., quantum systems, statis- 
tical systems, polymers, etc. (for a review see Ref. [pi, and 
references therein). 

DMRG is also a variational method and closely related 
to matrix product (MP) method In MP method 

the ground state of a system, which is assumed to be 
expressed as a matrix product, can be obtained by min- 
imizing the associated ground state energy density with 
respect to variational parameters. On the other hand 
DMRG method may find the ground state of a large sys- 
tem, which consists of Wilsonian blocks, B L • uB R , by 
keeping the most probable m states for each block B L I R . 
The most m probable states are selected by using the m- 
largest eigenvalues of the reduced density matrix p L l R 
for each block. For an overview and some connections 
to related fields including MP method, Ref. || is recom- 
mended. 

When a system of interest has a symmetry, the associ- 
ated eigenvalues of the density matrix p are degenerated 
and one should keep the all states that corresponding to 
the same eigenvalues. For example, if we consider a rota- 
tional invariant model, e.g., ID isotropic Heisenberg spin 
models, with the standard vertex-DMRG, we may use 
third components s z of spin as basis, then the eigenval- 
ues of the density matrix are degenerated due to the rota- 
tional symmetry. One thus needs to keep many states to 
improve numerical accuracy with s z basis in order to get 
a real physics at thermodynamic limit. There have been 
large-scale DMRG calculations with many states kept, 
for example, up to m = 300 @, m = 400 f§, m = 1700 
||, to estimate the 5 — 2 Haldane gap. However the 
more m states are kept, the more computational cost is 
necessary since the dimension of the Hilbert space for the 
superblock B L • uB R is proportional to m 2 . 



Sierra and Nishino (To) had applied the DMRG method 
to interaction round a face (IRF) models and developed 
the IRF-DMRG method. The advantage of using IRF- 
DMRG method in a system which has rotational sym- 
metry is that the dimensions of the associated Hilbert 
spaces are much smaller than that in standard vertex- 
DMRG, since the degeneracy due to the symmetry has 
been eliminated. They had demonstrated the power of 
the IRF-DMRG by calculating the ground state energies 
of the solid on solid (SOS) model, which is equivalent 
to spin-1/2 Heisenberg chain, and that of the restricted 
SOS (RSOS) model, which is equivalent to the quantum 
group invariant XX Z chain. They also had suggested a 
promising potential of the IRF-DMRG when it applies 
to higher integer quantum spin chains. Such a work has 
not yet been done as far as I know. 

In this work, the IRF-DMRG method is reviewed and 
further developed for the ID integer spin antiferromag- 
netic Heisenberg (AFH) quantum spin chains. In IRF 
formulation the dynamics of a model can be described 
by a local plaquette operator X$, which operates to the 
ith site of an IRF state as a "diagonal to diagonal trans- 
fer matrix" and their matrix elements are called the IRF 
weights as will be explained in the next section. We hence 
need the explicit expressions of the IRF weights for the 
higher integer quantum spin chain models in order to 
work with the IRF-DMRG. 

The rest of the paper is organized as follows. In the 
next section the IRF-DMRG method is reviewed. Af- 
ter the explanation of IRF formulation, the infinite sys- 
tem algorithm of the IRF-DMRG is discussed. How to 
target the excited states of AFH quantum spin chains 
are explained and the superblock configuration suited 
for targeting the excited state of the AFH spin chain 
is proposed. In Sec. Ill, the power of the IRF-DMRG is 
demonstrated by applying it to both 5=1 and 5 = 2 
AFH quantum spin chains. The Haldane gaps are cal- 
culated by keeping moderate number of m states. Fi- 
nally conclusions are stated. With the help of the Wigner 
Ekart's theorem summarized in Appendix A, I've derived 
the expression of the IRF weights for a nearest neighbor 
spin-5 interaction Si ■ Si+i in Appendix B. These expres- 
sions enable us to work with the IRF-DMRG. 
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FIG. 1. Diagrammatic representation of the operation of 
a plaquette operator Xi on a IRF state \£q, ■ • • ,£;, • • • ,£n), 
which is drawn as a zigzag chain of the lattice sites (open 
circles). The closed circle stands for the summation over the 
lattice variable ii. 



II. REVIEW OF IRF-DMRG 

The IRF-DMRG method is briefly reviewed here ac- 
cording to the Sierra and Nishino's paper ficf| . It is in- 
structive to begin with a brief review of graph IRF model 
in somewhat general context. The IRF representation 
of a rotational invariant quantum spin chain is then ex- 
plained with the help of the associated spin graphs. In 
fact ID quantum spin chain models can be formulated 
as a special case of graph IRF models. The algorithm of 
the IRF-DMRG is also reviewed. A method how to tar- 
get the excited states of AFH quantum spin chains and 
the appropriate superblock configurations are explained. 

A. Graph IRF model 

In interaction round a face or face models for short 
PlP^f , a state is represented with the lattice variables I 
assigned to the lattice sites; 

\£) = \£ ,£ u ...,e N ), (1) 

while the associated interaction is defined on a face, or 
plaquette of the lattice sites; 

Xi\ . . . , £i_i, l\, £ i+ i, . . .) 
= Vi?(£ J _ 1 1 t i+l \\...A-iAA+i,..): (2) 

where R are so called IRF weights, which play an impor- 
tant role such as Boltzmann weights in usual statistical 
models. The dynamics of the model is described by a 
plaquette operator Xi, which consists of the associated 
IRF weights. The plaquette operator Xi can be viewed 
as a local "diagonal transfer matrix" as shown in Fig. [j] 
in contrast with a usual "row to row transfer matrix." 

IRF model is characterized with a selection rule which 
determines whether the adjacent lattice sites of a lat- 
tice site are admissible or not. The selection rule can be 
described by its associated incident matrix Xg^>, whose 
elements are all either or 1. If A^/ = 0, then the 
lattice variables £ and £' cannot be assigned to adjacent 



lattice sites, i.e., not admissible to each other. The A 
characterizes the set of all possible configurations which 
contribute to the Hamiltonian or partition function of the 
IRF model. A configuration not derived from A neces- 
sarily has zero IRF weight. It may be more convenient to 
use a graph or diagram associated with the selection ride, 
or A, of an IRF model and this is called a graph IRF (or 
graph face) model. Each vertex of such a graph repre- 
sents a lattice site ii, which may take values from 1 to n,-, 
and they are connected with line if they are admissible. 

j 
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FIG. 2. Spin diagram for 5 = 1 chain. Each vertex repre- 
sents an IRF lattice site. The height of each site represents 
the total spin angular momentum ji summing from the most 
left vertex * (vacuum state) until the ith spin. Admissible 
vertices are connected with lines. 

For a rotational invariant spin chain, the correspond- 
ing selection rule is nothing but the addition rule of spin 
angular momenta, each spin state is hence classified with 
the total spin angular momentum j. The states of the 
chain with N spins will be given by the set \jw,M) 
(j mm < ]n < j max ) of states with total spin jx and the 
corresponding third component M . One of the graphi- 
cal representations associated with the addition rule of 
spin is the so called spin diagram. For example, the spin 
diagram for S = 1 spin chain is shown in Fig. ^. The 
corresponding lattice variables li of the spin diagram are 
the magnitudes ji of spin angular momenta until the ith 
spin. The height of each vertex represents the sum ji 
of the spin angular momenta starting from the most left 
spin until the ith spin. The ith and (i + l)-th vertices are 
connected if they are admissible, i.e., if they are satisfied 
with the addition rule; \ji — S\ < ji+i < ji + S, where 
S = 1 for spin-1 chain. Note that ji takes values from 
jf 1 ™ to jj max , i.e., the range of ji is site dependent. An 
IRF state of the N spin chain is represented by a path in 
the spin diagram |*, j\,ji, ■ • • ,ji, ■ ■ ■ , jjq) and the most 
left vertex * is corresponding to a vacuum state |0). If 
we know the all IRF weights associated with a spin chain 
model, we can work in the IRF formulation. It is hence 
important to derive the explicit expressions for the IRF 
weights of the spin chain model to apply the IRF-DMRG. 
In Appendix B, I have derived the expression of the IRF 
weights for the nearest neighbor spin-S* interaction and 
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summarized the diagrammatic representations and the 
corresponding IRF weights as a function of spin j for the 
5 = 1 isotropic AFH case in Fig. O. 



B. Algorithm 

The DMRG is an iterative algorithm to approximate a 
target state (the ground or excited state) of a large sys- 
tem which consists of the two blocks and lattice site(s). 
In each iteration each block is extended by adding one 
lattice site but keeping only the most m probable states, 
which are selected with the density matrix constructed 
from the target state. The superblock jl5[ is formed by 
the left block B L , the middle block which consists of a 
lattice point • or of two points ••, and the right block 
B R , which is not necessarily the reflection of B L . 



The superblock Hamiltonian H B * B can be con- 
structed from the two block Hamiltonians H B , H B 
and the IRF weights R using the following relation, 




(4) 



Figure |j shows the diagrammatic representation of the 
above equation. 

We then may find the ground (or excited) state of the 
superblock by using a Lanczos or similar method, 
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N = 2(k-l) + 4 



N = 2k + 4 



l*G>= ^ C, n j6.®&®»?c). (5) 

The left density matrix p BL * is readily constructed as 

" b ) = I>L.^U- ( 6 ) 
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FIG. 3. Schematic diagram of the infinite IRF-DMRG al- 
gorithm. Each block (dotted line box) is extended by adding 
single middle lattice point • in each DMRG iteration. Thick 
lines represent renormalized block states. 

Let me first review the infinite system algorithm |lC| ] 
of the IRF-DMRG. A Hilbert space of the superblock 
B L • B R can be written as 



The right density matrix p B * is constructed in the sim- 
ilar way if necessary. 

The next step is to diagonalize p B * to obtain the 
eigenvalues w b and the associated eigenvectors \u b ^J. 
Keeping the m states associated with m-largest w h (m = 
mj) in order to form the projection operator T'T. 



T = J2 Tb 



(7) 



H 



B mB 



= {\£a ® b®rj c 



A a ,b = Ab jC = 1}, (3) The operator T truncates the Hilbert space Ti B * into 



i.e., a space spanned by a left block state |£ ), a middle 
lattice point state |6), and a right block state \rj c ) (see 
Fig. ||), where the three lattice points (a, b, c) must be 
connected in a associated spin diagram. These lattice 
points take values from jf 1111 to j™ ax , where i = a,b, c. 
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FIG. 4. Diagrammatic representation for the construction 
the super block Hamiltonian 
responding to each term in Eq. 



of the super block Hamiltonian H Bl ' Br ; each diagram is cor- 



TL B , where B' L represents a block with one more lattice 
site than the block B L . A new left block Hamiltonian is 
then formed by using the projection operator as H B = 
T{H bL% )T^. In each IRF-DMRG iteration the system is 
extended by adding single middle site • to both blocks 
as shown in Fig. |[ 

The infinite system algorithm for the IRF-DMRG with 
the superblock B L • B R is summarized as follows: 

i) find the ground state of the superblock Hamilto- 
nian, Eq. (f|; 

ii) form the left (right) reduced density matrix p BL ' 
(p BR '), Eq. §; 

iii) diagonalize p BL * ( P bR ') to obtain the eigenvalues 
w h and the associated eigenvectors \ul) ( |ujj c )); 

iv) keep the m states |w| q ) corresponding to the m 
largest eigenvectors w b to form the projection op- 
erator T, Eq. (0); 
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v) renormalize the block operators by using T: 
H B ' L = T(H BL ')T\ etc.; 

vi) extend each block by adding one site •; 
repeat the processes i)-vi) for new blocks. 




FIG. 5. (a) Spin diagram of the superblock for targeting 
the excited state or the ground state of the S — 1 open spin 
chain ended with 5=1/2 spins, (b) The corresponding spin 
chain of N = 4 spin-1 spins (large hatched circles) ended with 
spin- 1/2 spins (small open circles). For targeting the excited 
states, which has total spin S T = 1, ferromagnetic coupled 
two spin-1/2 spins are attached. 



(a) 



(b) 




FIG. 6. (a) Spin diagram of S = 1 chain ended with 
S — 1/2 spins. AKLT state is corresponding to the hori- 
zontal line at j = 1/2. (b) The corresponding spin chain of 
N = 6 spin-1 spins (large hatched circles) ended with spin-1/2 
spins (small open circles). 
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FIG. 7. Schematic diagram of the proposed infinite 
IRF-DMRG algorithm. The left (right) block is extended by 
adding single middle left (right) point • in each IRF-DMRG 
iteration 



Next let me explain how to target the excited state 
of the AFH spin chain before discussing the superblock 
configuration suited for the excited state. Figure || shows 
how to target the excited state of the AFH spin-1 spin 
chain ended with spin-1/2 spins. In order to fix the to- 
tal spin St of the chain in St = 1 sector, two spin-1/2 
spins that ferromagnetically coupled with the coupling 
constant Jf < are attached as shown in Fig. g (b). 
Note that the attached two spins, which energetically fa- 
vor a triplet state, are not coupled with the chains. Since 
the whole system (the chain and the two spins) is set to a 
singlet state (0 sector), the triplet state of the two spins 
enforce St = 1 on the chain |l3| . 

Now let us consider the superblock configuration that 
suited for targeting the excited state. In the infinite sys- 
tem algorithm of the IRF-DMRG with the superblock 
configuration of B L • B R , the middle point is single 
bare site and plays an important role to successively im- 
prove the renormalized block Hamiltonians. In each IRF- 
DMRG iteration both blocks are renormalized with the 
middle site as shown in Fig. 0. In Ref. |l(]] the superblock 
configuration of B L • B R was used and it is suited for tar- 
geting the ground states of AFH spin chains, because the 
highest j site always lies in the middle of the superblock 
as shown in Fig. |q. Since the ground state of an AFH 
quantum spin chain with even N spins lies in the total 
spin sector, Jn = and obviously the vacuum state 
* has zero spin, i.e., jo = 0. Furthermore the system is 
reflection symmetric with respect to the middle point in 
the spin diagram. Thus the middle point always has the 
highest j in the IRF-DMRG formulation when targeting 
the ground state of the AFH spin chain. In other word 
the highest j of the left block, which is constructed from 
the left vacuum by adding spins, coincides with that of 
the right block, which is constructed from the right vac- 
uum, at the middle point in the spin graph. 




FIG. 8. Diagrammatic representation for the construction 
of the super block Hamiltonian H Bl " Br . 

However the situation is somewhat different for target- 
ing the excited state which lies in St = 1 sector. From 
Fig. [5] (a), we noticed that no single middle point at 
which the highest j of the left block coincides with that 
of the right one and that there is no reflection symmetry. 
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Hence it is better to use the superblock configuration 
of B L • »B R , at least when we use the infinite system 
method, and renormalize B L (B R ) with the middle left 
(right) site individually as shown in Fig. 0. The associ- 
ated superblock Hamiltonian can be constructed in the 
way diagrammatically shown in Fig. @. 



III. APPLICATION TO S = 1 AND S = 2 
ANTIFERROMAGNETIC HEISENBERG 
QUANTUM SPIN CHAINS 

The supremity of the IRF-DMRG method is demon- 
strated by applying to both 5=1 and 5 = 2 quantum 
spin chains. Haldane's conjecture jl4j that the physics 
of isotropic antiferromagnetic quantum spin chains de- 
pend substantially on whether the spin is integer or half- 
integer, has been motivating many physicists to study 
quantum spin chains. 5=1 AFH chain have been widely 
studied with various methods |16|. It is a well known 
fact that the ground state of an open 5=1 quantum 
spin chain has an effective 5 = 1/2 spin at each end. 
White |L7| had obtained the ground state energy per site 
of e = -1.401484038971(4) with to = 180 states kept 
and the Haldane gap of Ai = 0.41050(2) with to = 160 
states kept by using the standard s z -base DMRG. 




(b) 

FIG. 9. (a) Spin diagram of open S = 2 spin chain ended 
with S = 1 spins. The total S T = 1 state is targeted to find 
the excited state, (b) The corresponding spin chain of N = 4 
spin-2 spins (large closed circles) ended with spin-1 spins 
(hatched circles). The attached two spin-1/2 spins (small 
open circles) are coupled ferromagnetic or antiferromagnetic 
depending on the target state. 

5 = 2 AFH quantum spin chains have also been stud- 
ied 0j8|, although the numerical calculations are much 
elaborative than that of 5 = 1 case due to the longer 
correlation length £ and due to the larger number of de- 
grees of freedom per spin because the degeneracy due to 
the spin symmetry. The longer correlation length means 
that much longer chains (more than one thousand) are 
required to reach the convergence regime to exclude the 



finite size effects. Schollwock et al. estimated that 
the gap of A 2 = 0.085(5) and the ground state energy 
density to be eo = —4.761248(1) by using the standard 
DMRG with up to m — 210 states kept. Wang et al. 
H systematically analyzed and estimated that the gap 
A 2 = 0.0876 ± 0.0013 with up to m = 400 states kept 
during DMRG calculations. 

Following Schollwock et al, I also consider the open 
spin-5 AFH quantum spin chain terminated with spin- 
5/2 spins to cancel out the effective spin-5/2 spins: 



H — JendSl • S 2 



N-2 

J 2_j ^ ' + JendSjV-l ' Sat, (8) 



where the both coupling constants J and J en d are set to 
unity for simplicity. The spin diagram for the 5 = 1 
ground state is shown in Fig. |(| I have obtained the 
comparable result of e = -1.40148403897 with the IRF- 
DMRG by keeping only m — 80 states, which consist 
of 25(j = 1/2), 31(j = 3/2), 19(j = 5/2), and 5(j = 
7/2) rrij states, without resorting to a scaling technique. 
Note that since rrij states in j base correspond to (2j + 
l)m,j in s z base, the above to = 80 correspond to to = 
328 in s z base. For the 5 = 2 chain with N w 200, I 
have found that the ground state energy density is eo = 
—4.7612481(6) by keeping to = 90 states, which consist 
of 9(j = 0), 24(j = 1), 27(j = 2), 19(j = 3), 9(j = 4), 
and 2(j = 5) rrij states, hence corresponding to to = 452 
states in s z base. 

The finite size correction to the Haldane gap As with 
open boundary conditions is proportional to the inverse 
of the square chain length according to the ID field the- 
ory 



A(m = oo, N) = Ac 



2^ + °<^>- < 9 > 



where v is the spin wave velocity and As the Haldane 
gap of the spin-5 AFH spin chain at the thermodynamic 
limit. Figure [Tc| indicates the gap Ai(m, TV) as measured 
by the difference between the lowest energy of 5 T = 1 
states and that of 5 T = states, as a function of the 
spin chain length N and the number of states to kept 
in the IRF-DMRG iterations. The excited states were 
obtained by using the method explained in the previous 
section as shown in Figs. || and ^. With the simultaneous 
extrapolation for to and 1/TV, the Haldane gap for the 
5 = 1 AFH chain is estimated to be Ai ^ 0.4104(5). 

Wang et al. [|| pointed out that one cannot use the ex- 
trapolation with respect to 1/N to obtain the gap value 
at the thermodynamic limit when to is not sufficiently 
large: the minimum of each curve for As(m, N) deviates 
from the vertical axis as to decreases; while Eq. (^) tells 
us the minimum should be located just on the vertical 
axis in the limit to — > oo. DMRG involves a systematic 
error associated with the truncation or keeping a finite 
to states of the renormalized Hilbert space. They thus 
remarked that these errors are not so serious for 5 = 1/2 



5 



or S = 1 AFH chains but the errors become crucial for 
higher spin chains and scaling for m should be carefully 
carried out. In standard DMRG calculations the multi- 
plicity 2S + 1 of spin-5 is not eliminated. Then in order 
to treat the higher spin chain, the larger number m of 
states should be kept. The IRF-DMRG however has a 
great merit to eliminate the degeneracy due to the spin 
symmetry. We can hence keep effectively larger m states! 
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FIG. 10. The gap A 1 (m.,iV) in the unit of the coupling 
constant J as a function of the number of state m = ^ . nij 
kept in the IRF-DMRG calculations and the spin chain length 
N. Each cross denotes the position of minimum for a given 
m curve. 



and S = 2 isotropic AFH quantum spin chains. With 
the moderate number (up to m = 90) of states kept in 
the IRF-DMRG iterations, the Haldane gaps and ground 
state energy densities were readily calculated since the 
degeneracy due to the spin symmetry had been elimi- 
nated. 
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FIG. 11. The gap A 2 (m,N) in the unit of the coupling 
constant J as a function of the number of state m = y\ rrij 
kept in the IRF-DMRG calculations and the spin chain length 
N. Each cross denotes the position of minimum for a given m 
curve. The open and solid diamonds denote upper estimated 
value A(j and lower one A_l, respectively. 



Figure [ll] shows the gap A2 (m,N), as measured by 
the difference between the lowest energy of S T = 1 states 
and that of S T = states, as a function of the spin chain 
length N and the number of states m kept in the IRF- 
DMRG iterations. The extrapolations are performed for 
both m and N simultaneously and the upper estimated 
value Ajj and lower one Al are obtained by using the 
polynomial fits with second order and third order, re- 
spectively. The estimated Haldane gap for the S = 2 
AFH chain is A 2 ^ 0.0878 ± 0.0016. 



IV. CONCLUSIONS 

The IRF-DMRG is reviewed and developed for higher 
integer quantum spin chain models which has a rotational 
symmetry. The explicit expressions of the IRF weights 
for the nearest neighbor spin-S 1 interaction have been 
derived as a function of total spin j. Using these IRF 
weights the IRF-DMRG has been applied to both 5=1 
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APPENDIX A: WIGNER-ECKART'S THEOREM 

The Wigner-Eckart's theorem is briefly reviewed to de- 
rive the IRF weights associated with general spin-S* in- 
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teraction in the form of S • S' in the next appendix. Most (SMlT^lSM 1 ) — (— l) s ~ RI SIS 
of the results presented here are already known and can —MO M 

be found in Ref. M. 

(k) 

Let T^j is an irreducible tensor operator with an 
angular momentum k and the third component M(= Since the left hand side is MSm m>, one gets 

.(*0 



— fc, • • • , fc). Tjrf commutes with the total angular mo- 
mentum operator J of a system considered as 



(S||S||5) = y/S(S+l)(2S+l). 



(S\\S\\S). 
(A8) 

(A9) 



[Jz,T 



M J 



MT 



(fc) 



(Al) 



[J x ± z J y , T$] = y/k(k + l)-M(M±l)T$ ±1 . (A2) 



The following relations Q are used to derive the ex- 
pression for the IRF-weights of general spin-5 1 chain in 
the next appendix: 

W T,(fe)|„./ 



i.e., is transformed as a tensor under a rotational (jihJM\T 1 ■ T 2 \j' x j' 2 J'M 



i.e., J-^ is u tuitsiui muu as a lciisui unuci a louauunai 

operation for the system. For example, each component 
of spin operator S is expressed as 



T (i) _ o 



(A3) 



The inner product of a pair of irreducible tensors T^ k ' 
and U^ k ' is defined with 



T (k).u<k)= £ 

M=-k 



-Mrp{h) TT {k) 
A M u —M' 



(A4) 



For spin operator one readily checks the following rela- 
tion: 



T (i) . T (i) _ g . 



(A5) 



<W«W(-i) i # # fc 



x(ii||Ti fc) ||iO(i2||Tr J ||^) 



(fc) 1 1 " 



(j^JIIT^Hj^J') 
= W" 1 



ji+h+J'+k J jl J J2 

J' # fc 



(iii2J||Tf||ibV) 



< v /(2J+l)(2J' + l)0 2 ||T«||^; 



(A10) 



xV(2J+l)(2J' + l)( J1 ||TW|| Jl ), (All) 



(A12) 



Wigner-Eckart's theorem shows the matrix elements 
for are factored as a product of a configuration de- 
pendent part, which expressed as Wigner's 3-j symbols 
or Clebsh-Gordan (CG) coefficients, and a configuration 
independent part: 



(JM\T {k) \j'M r 



= (-1 



J k J' 
-M /i M' 



(J||T«||J'), 



(A6) 



where (J\\T^\\J') are called reduced matrix elements, 
which is independent of the configuration M or M' . 
Wigner's 3-j symbol is related with CG coefficient as 



J k J' 
-M n M 



(-1 



J-k-M' 



V2J 7 



(J —Mkfj,\J' —M') 



(A7) 



Once the reduced matrix elements are known, one gets 
the all matrix elements with only calculations for CG 
coefficients. For example one easily finds the reduced 
matrix ele ment for spin operator S by applying the the- 
orem Eq. ( |A6| ) to S z , which is corresponding to fc = 1 
and /i = 0, as the following: 



APPENDIX B: IRF-WEIGHTS 

I here derive the IRF weights for the following Hamil- 
tonian H, which is invariant under rotations. 



Hi,i+1 — Sj • Sj_| 



(Bl) 



where Si and Sj+i are not necessarily same. The IRF 
weights for the H are expressed as the matrix elements 

((JiSi+i), J i+1 M i+1 \Si ■ S i+1 \(4S i+1 ), J' l+1 M' l+1 ), (B2) 

where |(«7jiSj+i), Jj+i-Mj+i) is a state that the sum of 
spin angular momenta until the (i + l)-th spin is J;+i 
and the corresponding third component is Mi+i- The 
J;+i consists of Jj and a spin-Si+i spin. 
Using Eq. (A10) one finds 



{{JiSi+i), Ji+\Mi + i\Si ■ Si+iKJ^Sj+i), Ji_|_iM,+i) 

\j i+1 +Ji+s 1+1 | Ji S i+ i J l+ \ 
\ S i+i J[ 1 

x(S i+ i||S||S i+ i)(.7 i ||S||.7!). (B3) 



= (-I)" 



Since |Jj) can be expressed as a tensor product of Ji-\ 
and spin-S*, the last reduced matrix element in t he ab ove 
equation may be written, with the help of Eq. (All), as 
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(j,||s||j;) = Vi-xSiJiWawjuSiJl) 

_f , f -i\Ji+Si+Ji-i+l j Si Ji Ji-l 

-°J<-iJU( L ) 1 J' Si 1 



x v /(2J i + l)(2J? + l)(5 i ||S||5 i ). (B4) 
Hence the final expression is obtained as 
J- 



R I Ji-i Ji+i 
J,. 



— ((JiSi+i), Jj+i|Sj • Sj+i|(J i S'i+i), Ji+i) 



^ySi+S^+J^+Ji+Jl+J^+l h 2J . + x ^j, + ^ 



Xy/Si(Si + l)(2S t + l)^/S~^XSi^ 1 + l)(2S i+1 + 1) 

Ji Sj+x Ji+l 1 J Si Ji Ji-l 

Si+i Ja 1 j 1 Ji S*i 1 
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(B5) [6 

[7; 



1. S = 1 nearest-neighbor interaction 

There are thirteen nontrivial IRF weights R for the 
5 = 1 nearest neighbor spin interactions: Sj • They 
can be obtained by substituting Si = = 1 in Eq. 

(B5) and the results as a function of j are summarized 
in Fig. [l2] with diagrammatic representation. In each di- 
agram the magnitude of spin angular momentum for the 
middle vertex is assumed to be j, and the height of each 
vertex represents the magnitude of the spin angular mo- 
mentum that assigned to the vertex. The hrst diagram, 

( 3 \ 

for example, stands for R I j + 1 j + 1 I = -q^-. 

V 3 J 

Similar results for 5 = 2 case are readily obtained with 

the help of a symbolic manipulation language such as 

MATHEMATICA. 



V , + i 


1 

V. 7 + 1 


\ 1 




~ 7(7+1) 


-1 


/ 1 




.7 





.7+2 /2.7+1 




ViO+2) 


J + 1V2J+3 


i+i 


J7 


Vi(i+2) 




.7 /2.J+3 


i+i 


O'+DV 2 ^ 1 



FIG. 12. The diagrammatic representations and corre- 
sponding expressions for the IRF-weights of the S = 1 near- 
est-neighbor interaction Si ■ Si+i. In each diagram the hight 
of a vertex represents the magnitude of the spin angular mo- 
mentum assigned to the vertex. The height of the (lower) 
middle vertex is assumed to be j. 
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